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^ : ABSTRACT 

2 I Measurements of disoriented chiral condensates in heavy ion collisions 

Q\ • at RHIC and the LHC can yield fundamental information on the nature of 

fH ■ the QCD phase transition. I review theoretical efforts to understand the 

Qh! evolution of the condensate and present new results on experimental signals 

p ' in the single pion spectrum and in pion interferometry. 



Invited talk at the Third International Workshop on 

Relativistic Aspects of Nuclear Physics, 

24-27 August 1993, Rio de Janeiro, Brazil. 



This manuscript has been authored under contract number DE-AC02-76CH00016 with the U.S. Depart- 
ment of Energy. Accordingly, the U.S. Government retains a non-exclusive, royalty-free license to publish or 
reproduce the published form of this contribution, or allow others to do so, for U.S. Government purposes. 



The order of the chiral restoration phase transition in high temperature QCD is cur- 
rently unknown. If the transition is second order, Rajagopal and Wilczek (RW) propose 
that the nonequihbrium dynamics in ion-ion coUisions at RHIC can generate transient 
domains in which a macroscopic pion field develops.^ These domains are called disoriented 
chiral condensates (DCC). Bjorken, Kowalski, Taylor and others pointed out that DCCs 
can lead to fluctuations in the charged and neutral pion spectra.^ Hints of this behavior 
may have been seen in the Centauro cosmic ray events.^ In ion-ion collisions, the ability of 
experimenters to identify DCCs amidst the background produced by conventional particle 
production and scattering critically depends on the domains' size and energy content. 

In this talk, I explore the RW mechanism for domain formation to determine the size 
of domains in nuclear collisions, reporting on work with Gocksch and Pisarski^ and Miiller 
^. I then briefly speculate on phenomenological consequences of DCC formation. [For a 
discussion of J/'^ suppression as presented at this meeting, see Ref. [6].] 

Equilibrium high temperature QCD manifests a chiral symmetry if the light up and 
down quarks are taken to be massless. However, a phase transition occurs at a critical 
temperature T^ ~ 140 MeV in which chiral symmetry is broken by the formation of a scalar 
(iqq) condensate. In the real world chiral symmetry is explicitly broken at all temperatures 
by the few-MeV current quark masses. In that case the nature of the transition has yet 
to be established.^ [For reviews of lattice and flnite temperature QCD, see e.g. Refs. [8].] 

Rajagopal and Wilczek pointed out that the chiral condensate can become temporarily 
disoriented in the nonequihbrium conditions encountered in heavy ion collisions. Near Tg, 
the approximate chiral symmetry implies that the scalar condensate is nearly equivalent 
to a pion-like pseudoscalar isovector condensate ~ (^75^^) , where r are the Pauli isospin 
matrices. Consequently, domains containing a macroscopic pion fleld can appear as the 
temperature drops below T^. Such domains will eventually disappear as the system evolves 
towards the true vacuum in which only the scalar condensate is nonzero. In the flnite-sized 
heavy ion system, however, the evolving domains can radiate pions preferentially according 
to their isospin content. 

The excitation of a pion fleld over a substantial fraction of the collision volume can 
result in novel fluctuations in the number of neutral and charged pions. Consider a single 
ideal domain in which the pion fleld is oscillating along a flxed isospin direction yf. Pion 
production from this domain is proportional to the square of the fleld strength, so that the 
fraction of neutral pions / is 

/ = KT/ E 1-1' - ^°^' ^- (1) 

0,+,- 

AU orientations of tt are equally likely, so that the probability of a flnding a fraction / is 

P(/)ocrf(cos e)/df^r^'^ (2) 

It is therefore most probable for a single domain to emit more charged pions than neutral.^ 
The detectability of domains in a nuclear collision depends on their size.^ If the inter- 
action volume in a nuclear collision is dominated by a single domain, then DCC formation 
can lead to a measurable isospin asymmetry. If a nuclear collision instead produces many 



uncorrelated small domains, the spectrum of fluctuations would be gaussian, as expected, 
e.g., if the pions were produced independently. 

To develop a more concrete picture of how the chiral condensate evolves, we use the 
linear sigma model, in which the pion field is coupled to a scalar a field that characterizes 
the scalar condensate . The fields interact through a potential 

V = \[^^ + a'^-v^f /A-Ha. (3) 

Chiral symmetry transformations correspond to 0(4) rotations of the vector field $ = 
((7,7?). For H = this symmetry is spontaneously broken at zero temperature, with 
{a) — V. Current algebra implies that {a) = /jr = 93 MeV.^ The external field H (oc the 
light quark masses) breaks the 0(4) symmetry explicity and gives the pions a small mass 
TOtj- = ^JH/ Jt^ = 140 MeV. In comparison, the sigma mass 171^ ~ v2\v^ is quite large, 
perhaps ~ 600 MeV. 

Following Ref. [1], I will take this model with the parameters fixed at T = as the 
basis for a Ginzburg-Landau description of the chiral transition. For H = 0, the transition 
would then be second order with an 0(4) order parameter $. DCC domains would then 
be well defined entities in the thermodynamic limit. While there is strictly no transition 
for H ^ 0, transient domains can form as follows. 

In the idealized 'quench' scenario proposed by Rajagopal and Wilczek, one assumes 
that the system initially has $ = 0, as appropriate at high temperature, and follows the 
development of the system using T = equations of motion derived from V. The initial 
state of the system is clearly unstable. The system "rolls down" from the unstable local 
maximum of V{^) towards the nearly stable values with |$| = t; (the symmetry breaking 
term —Ha is relatively small). Field configurations with n ^ develop during the roll- 
down period. The field will eventually settle into stable oscillations about the unique 
vacuum (/jr, 0) for if 7^ 0, but oscillations continue until interactions eventually damp the 
motion. 

To be more concrete, the linearized equations of motion for the Fourier components of 
the pion field are: 

i^ffj={A„^-An*j. (4) 

Field configurations with ($) = and momentum k < vAw are unstable and grow expo- 
nentially; modes with higher momenta do not grow. The k = mode grows the fastest, 
with a time scale 

TR = {At;2}-i/2 ^ ^/^^. (5) 

The unstable modes grow for a time of order tji until $^ approaches v'^, after which the 
system oscillates about a = f-^-. 

Significantly, the time during which growth can occur is quite short, tr ~ 0.5 fni for 
rria ~ 600 MeV. Rajagopal and Wilczek found that the power oc (tt^^tt^^^) in the low 
momentum pion modes grows when the exact classical equations of motion are integrated, 
demonstrating that domains indeed form. [Here (. . .) represents an average over the initial 
of the fields, which are taken to have a thermal spectrum.] What is not clear, however, is 



whether the time during which the system is unstable is sufficient for truly large domains 
to form. 

Gocksch, Pisarski and I^ studied the domain size by numerically integrating the equa- 
tions of motion in the quench scenario and extracting the spatial correlation function: 

(tt (f , t) n (0, t)) oc J d3A;(7r_.7r^)e-^^-^ (6) 

where isospin labels are implicit. This correlation function contains size information since 
(7r(x, t)7r(0, t)) ~ (7r(0, t))^ holds inside the domain. The numerical results in Ref. [4] 
can be understood by taking TTk{t) from the linearized equation (4) and assuming that 
the fluctuations are initially thermal, so that (tt rTrr) ~ {E{exp E/Tc — 1)}^^ for E = 

{k +m^} ' . A saddle-point integration yields (7r(x)7r(0)) oc e^^ /s^fl^. The domain size, 
Rd ~ 2ri? oc m~^, is set by the sigma mass. For ma = 600 MeV, we find that domains are 
essentially pion sized. Numerical simulations confirm this result. 

The lesson drawn from the quench scenario is that large nucleus-sized domains can 
occur only if the sigma mass is small. ^'^'^'^^ However, it is quite possible that the effective 
nicr is reduced in the high energy density heavy ion environment. In fact, at a truly second 
order phase transition we expect rricr — the inverse correlation length — to vanish! In that 
case, the parameters of the effective potential for our Ginzburg-Landau model cannot be 
naively taken from zero-temperature physics as we have done. 

To explore the role of the medium in domain formation, Miiller and I studied the 
evolution of the condensate in the presence of a nonequilibrium bath of quasiparticles.^ 
We find that larger domains are possible due to two effects, i) the reduction of ma and ii) 
the "annealing" of the system, rather than quenching, by the slowly expanding heavy ion 
system. [Our use of the word "annealing" is meant to distinguish our scenario from RW's 
"quench" — the similarity of these terms to those used in condensed matter physics is in 
some respects misleading.] 

I illustrate the role of mass reduction by considering a system near equilibrium. The 
system is described by an effective potential V^ff that has the behavior shown in Fig. 1.^^ In 
the Hartree mean-field approximation, the change of Ves{^) as a function of temperature 
is determined essentially by the coefficient of $^ , i.e. the mass term. Taking i7 = 0, one 
obtains a linearized equation analogous to (4) that implies a growth rate r^ ^^ oc m^ oc 

(T^ — T^)^'^. Growth is very slow near the Tg ^ \/2v, prolonging the time over which 
the system is unstable. Only as the potential approaches its free space shape does the 
roll-down become rapid. 

The annealing scenario describes the nonequilibrium evolution of the system in a regime 
where the time scale Tev for the evolution of Vqs is much longer than tr. One expects 
such a slow cooling scenario to further enhance the domain size — in metallurgy, the 
more rapidly one quenches, the smaller are the resulting crystals. The evolution of V^efj is 
determined by the expansion of the surrounding quasiparticles, which constitute a time- 
varying nonequilibrium heat bath. The annealing regime is relevant at RHIC energies 
because the temperature falls to Tg ^ 140 MeV only at very late times, perhaps Tc ~ 
10 — 20 fm. The expansion time Tev is then '^ Tc^ tji^ 0.5 fm. 
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Figure 1: The effective potential for various temperatures. 

Miiller and I explored the evolution in the annealing regime using a nonequihbrium 
Hartree-like approximation. [See Ref. [5] for details.] We indeed find that anneahng 
prolongs the time during which the system is unstable compared to a quench. Domain 
sizes up to 7 fm are possible if we assume that m^{Tc) = {i.e., a second order transition). 
More conservatively, we find sizes of 3 — 4 fm for m^{Tc) = 300 MeV. Such large masses 
are possible, e.g., in the large A^ limit of the 0{N) model studied by Boyanovsky et a/.^^ 
and Kluger et al.^^ However, QCD may be far richer than the large A^ limit implies. 

Suppose that a nuclear collision produces a large domain, either in a typical event due 
to the annealing dynamics or due to a rare fiuctuation. How can we tell? The number of 
pions is not large, even ideally. The amount of energy available for pion production per 
unit volume is the potential difference AV^ ~ 60 MeV fm~"^ between $ = and $ = (/^r, 0). 
If all of the energy in a domain of size Rd ~ 3 — 6 fm goes into pion production, one expects 
A'dcc ~ AVR^/rrin ~ 30 — 300 pions. In comparison, conventional production mechanisms 
in a central Au+Au collision at RHIC can produce 1000 pions per unit rapidity; this 
constitutes a background to the DCC signal. 

The source of pions from this domain is coherent and, therefore, concentrated at low 
transverse momenta. Horn and Silver, Gyulassy et ai, and Amado et al. have developed 
a coherent state formalism applicable to DCC production. ^^ The amplitude for emitting 
a single neutral pion from a space time point (xj, tj) within a domain is real and satisfies 
Aj oc Ti^{xj) oc exp{— Xy/2i?|j}. The production rate from a single domain is 
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y^ Aje^^P'''3+^Ep^i\'^ ~ e 



^p'Rl 



(7) 



neglecting an energy-dependent factor that is slowly varying for small momenta. The rate 
is largest for p^ < R^ . 



The DCC signal (6) appears atop an isospin symmetric background. Assuming that aU 
the DCC pions are neutral and taking the background produced by the lund/fritiof event 
generator for central Au+Au, I obtain the spectra in Fig. 2. A signal that is significant 
compared to statistical fiuctuations in the background emerges at low pt for a domain of 
size Rd > 3. 
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Figure 2: Neutral pion production from a single large DCC in Au+Au 
at RHIC. 



Experimentally, one must search for DCCs on an event-by-event basis, looking for 
significant fiuctuations in the neutral and charged pions. A domain producing a neutral 
fraction / introduces a contribution fE{dN jSp^^^c to the tx^ spectrum, with a correspond- 
ing contribution {(1 — f)/2]E{dN /Sp)^^.^. to the tt"*" and 7r~ spectra. The distribution of 
/ in a sample of large-domain events is given by (1). 

To check that the pion emission from a candidate large-domain event is coherent, one 
can study the Bose correlations of pairs of identical pions. Pions produced by conventional 
scattering mechanisms are largely incoherent, and exhibit an intensity interference that 
is analogous to the Hanbury-Brown-Twiss effect for photons. ^^ Measurements of tx-k cor- 
relations are typically used to study information on the size of the interaction volume in 
nuclear collisions, much as the HBT effect is used to measure the size of stars (although 
interpretation of experiments is highly nontrivial^^). In contrast, a fully coherent source 
of pions exhibits interference at the amplitude level. Consequently, the single particle 
spectrum (7) from an isolated domain depends on the size of the domain, while the pair 
distribution is simply the product of single particle spectra. 

We therefore expect the HBT effect to be suppressed in a nuclear collision if a large 
domain forms. However, the HBT effect does not completely disappear, since intensity 



interference can occur when a pion from the DCC has a similar momentum to a background 
pion. To estimate this effect, I follow Gyulassy et al. in Ref. [13] and write the two pion 
correlation function 



C (pi,p2) = l + {l-D (pi)) {l-D {p2)) ~p {qf 

+ 2{D (pi) D {p2) {l-D (pi)) {l-D {P2))V^^P {q) 



(8) 



where q = pi — pi is the relative momentum of the pair and p is the Fourier transform of 
the space-time density of the background pions. Observe that (8) depends strongly on the 
fraction of coherent pions D = {dN / d^p)^^^^./ {dN / d?p)tot- 

The possible impact of a DCC on tt^tt^ correlations is illustrated in Fig. 3 as a function 
of domain size. I assume that the domain produces only charged pions and that the 
transverse source size for the fritiof background is 7 fni. A combination of one- and two- 
particle measurements can allow one to disentangle the DCC signal from the background. 
Of course, the interpretation of experiments will be extremely tricky! 
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Figure 3: tt^tt— correlations as a function of qT = Pti —PT2 from a single 
large DCC plus a background of incoherent pions in Au+Au at RHIC. The 
total momentum of the pair is taken to be K = {p^i +pr2)/2 = 0. 

The schematic results in Figs. 2 and 3 indicate that domains larger than 3 fm can have 
measurable consequences. If seen in experiments, what would they tell us about QCD? All 
of these signals — the isospin fluctuations, the enhancement of the pion spectrum at low 
Pt, and the suppression of HBT correlations — are characteristics of any large coherent 
source. In principle, Bose condensation in an ideal pion gas can produce similar signals. ^^ 
While it is argued^^ that the conditions at RHIC are not favorable for Bose condensation 
in the absence of a potential such as (3), systematic experimentation will nevertheless be 
needed to prove that the coherence comes from a disorientation of the chiral condensate. 
If an annealing scenario is valid, I expect large domains to occur in most central collisions. 



Lego plots^ could then reveal a 'clumpy' event structure indicative of different domains. 
The pion excess within each domain would be independent of the background multiplicity. 
On the other hand, the number of condensed pions in ideal Bose condensation would grow 
with the overall multiplicity and have a homogeneous structure. 

I am grateful to Andreas Gocksch, Rob Pisarski, and Berndt Mueller for their very 
enjoyable collaboration. 
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